Empirical Proton-Neutron Interactions and Nuclear Density Functional 
Theory: Global, Regional and Local Comparisons 
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Calculations of nuclear masses, using nuclear density functional theory, are presented for even- 
even nuclei spanning the nuclear chart. The resulting binding energy differences can be interpreted 
in terms of valence proton-neutron interactions. These are compared globally, regionally, and locally 
with empirical values. Overall, excellent agreement is obtained. Discrepancies highlight neglected 
degrees of freedom and can point to improved density functionals. 
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As with other many-body systems, the structure of the 
atomic nucleus depends on the interactions of its con- 
stituents, protons and neutrons. These interactions, re- 
flecting the strong and Coulomb forces, and the Pauli 
Principle, are complex. Nevertheless, their understand- 
ing is critical to interpreting nuclear structure and its 
evolution with N and Z . Similar issues arise in other fi- 
nite complex systems, such as nanostructures, and there 
is increasing overlap in the theoretical tools applied, fn 
nuclei, where two kinds of fermions come into play, the 
proton-neutron (p-n) interaction plays the key role in the 
development of long-range collective correlations, includ- 
ing non-spherical shapes. Due to the shell structure of 
nuclei, p-n interactions of the valence (open shell) nucle- 
ons are the most important. 

Since nuclear masses embody the sum of all nucleonic 
interactions, they provide a laboratory in which it is pos- 
sible to isolate and extract specific interactions using ap- 
propriate mass differences In particular, the average 
interaction of the last two protons with the last two neu- 
trons in an even-even nucleus is given by the following 
double difference of binding energies 0, Q : 



6V^n{Z,N)^ \[{B{Z,N)-B{Z,N-2)} 



{B[Z ~2,N)- B{Z -2,N - 2)}] 
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With the 2003 mass evaluation [ij , it became possible to 
evaluate a much larger set of 5Vpn values. These have re- 
vealed 0, d, 0, [1] striking bifurcations near closed shells 
and systematic patterns spanning major shells a cor- 
relation between SVpn values and growth rates of collec- 
tivity [sl; and intriguing patterns in specific regions 
While simple calculations with schematic zero-range in- 
teractions give reasonable results in the deformed rare 
earth nuclei, they fail completely in the actinides 
Clearly, a more sophisticated approach is needed. 



The indicator ([T]) involves masses of four neighboring 
even-even nuclei. Theoretical understanding of the be- 
havior of SVpn throughout the whole nuclear chart thus 
calls for an approach that is capable of predicting nu- 
clear masses with arbitrary Z, N values. Such an ap- 
proach must fulfill several strict requirements. First, 
it should be rooted in microscopic theory. Second, it 
must be general enough to be confidently applied to 
regions of the nuclear landscape whose properties are 
largely unknown. Third, it should be capable of han- 
dling symmetry-breaking effects resulting in a variety of 
intrinsic nuclear deformations. These requirements are 
met by density functional theory (DFT) in the formu- 
lation of Kohn and Sham 0. The main ingredient of 
the nuclear DFT [13, El is the energy density functional 
describing conditions locally around each nucleon. This 
can be realized by expressing the functional in terms of 
local nucleonic densities and currents. The energy func- 
tional is augmented by the pairing term describing nu- 
clear superfluidity fl^ . When not corrected by additional 
phenomenological terms, standard functionals, treated 
self-consistently, reproduce total bindin g en ergies with 
an rms error of 1.5 to 4 MeV 0, [ll, 111. However, 
they have been successfully tested over the whole nu- 
clear chart for a broad range of phenomena, and usually 
perform better when applied to energy differences and 
other global nuclear properties such as radii and nuclear 
deformations. 

As pointed out in Ref. TB^ , from Eq. ([T|) , 5Vpn approx- 
imates the mixed partial derivative 



SVpn{Z,N) 
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For nuclei with an appreciable neutron excess (Tj>l), the 
average value SVpn probes the symmetry energy term in 
the macroscopic mass formula 16, 17]: 



SV, 



-1/3 
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where agym and assym are the symmetry and surface- 
symmetry energy coefRcients, respectively. For nuclei 
with N ^ Z, SVpn also contains the Wigner energy which 
can be extracted by taking differences between values of 
SVpn in neighboring nuclei [3] ■ On top of SVpn , detailed 
fluctuations of SVpn carry important information about 
shell effects and many-body correlations. 

It is the purpose of this Letter to present the first 
results of large-scale microscopic calculations for SVpn 
within the DFT framework. Our work demonstrates 
that, since SVpn is a relative quantity, calculations with 
realistic functionals can reproduce empirical values ([T]) to 
significantly better than 100 keV in many mass regions. 
This is of the same magnitude as local fluctuations in em- 
pirical SVpn values. Hence, these calculations can assess 
the adequacy of currently used energy density functionals 
and help search for improvements in such aspects as their 
density dependence and dynamical correlation effects. 

The present large-scale calculations of nuclear masses 
are based on the HFB-hTHO code described in Ref. [H. 
Our DFT approach is based on a self-consistent solution 
of the Hartree-Fock-Bogoliubov equations with an ap- 
proximate Lipkin-Nogami treatment of pairing followed 
by an exact particle number projection. Calculations are 
performed in a transformed harmonic basis spanning 20 
major shells. In the particle-hole channel we employed 
the SkP [3 and SLy4 Skyrme functionals. The 
pairing functional corresponds to a density-dependent S 
interaction of Ref. [2l[. Unless otherwise indicated, a 
mixed-type pairing was employed. Overall, over 1000 
nuclei were studied. As SVpn corresponds to a second 
derivative, in order to obtain reliable values, calculations 
required very high numerical accuracy. In fact, we found 
out that SVpn provides an excellent check on the precision 
of calculations. (A similar conclusion was drawn earlier 
[§| in the context of the accuracy of measured masses.) 

We first present a global comparison of experimental 
and theoretical SVpn values in Fig. 1. To accommodate 
the wide range of SVpn values, while preserving the visi- 
bility of the microstructure for heavier nuclei, the panels 
use different vertical scales. The overall trends in the 
data are well reproduced by Eq. ([3]) with the SkP val- 
ues of asyin=30MeV and assyni=-45 MeV [13]. This con- 
firms empirically the importance of the surface symmetry 
term. Superposed on this secular decrease there are con- 
siderable fluctuations. The most dramatic effect is the 
empirical singularities in light N—Z nuclei. This is well 
understood [sl, as resulting in part from r=0 (p-n) 
pairing interactions which are not considered in our DFT 
model. The regional fluctuations on the right seem to 
have substantial differences between theory and experi- 
ment, but a more detailed view provides a more accurate 
perspective on regions of agreement and regions where 
the calculations are missing key ingredients. 

To pursue this. Fig. 2 shows four isotope chains com- 
prising vibrational systems (Cd), two deformed chains 
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FIG. 1: Comparison of empirical SVpn values (top) across the 
nuclear chart with DFT calculations using the SkP functional 
(bottom). The plot includes stable nuclei as well as available 
results for the trans-Pb nuclei from Po to U. The average SVpn 
values from Eq. (|3]) are shown by a solid gray curve. 

(Ra, U), and proton magic Pb nuclei. The calculations 
were carried out with the SkP functional with the stan- 
dard mixed pairing term (SkP) or with the volume pair- 
ing (SkPV). For Cd and U, the agreement with experi- 
ment is remarkably good for both functionals, and quite 
acceptable for Ra, generally within a few to lO's of keV. 
Given the average mass accuracy of several MeV, this 
demonstrates a striking ability to probe specific interac- 
tions by exploiting the filtering capabilities of mass differ- 
ence indicators. The Ra-U panels show opposite empiri- 
cal trends, both well reproduced by the calculations. Yet, 
there are discrepancies in Ra for iV= 128-134 which point 
to missing physics. Here, an octupole, reflection asym- 
metric, degree of freedom, not included in the present 
calculations, plays an important role (2^. The Pb nuclei 
show large discrepancies, reflecting the inadequacy of the 
current DFT approach for describing strong dynamical 
changes in "core" structure near magic numbers. 

Another interesting feature seen in the SkP and SkPV 
results of Fig. 2 is that the choice of pairing only weakly 
affects SVpn- This is because, at least for nuclei away from 
the N=Z line, the pairing correlation energy can roughly 
be written as a sum of independent proton (p-p) and neu- 
tron (n-n) contributions; hence, the leading components 
of pairing are expected to be filtered out by the indicator 
It is instructive to inspect the effectiveness of this 
filter. Figure 2e shows the differences for the SkP and 
SkPV calculations in masses, and SVpn- Although 
the mass differences {BE) range from up to nearly an 
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FIG. 2: Panels (a-d) display empirical and calculated SVp„ 
values for Cd, Pb, Ra, and U isotopes. Calculations are based 
on the SkP functional with mixed (SkP) or volume (SkPV) 
pairing. Panel (e) shows differences in binding energies (BE), 
2-nucleon separation energies {S2n) and SVpn values obtained 
in SkPV and SkP for the Cd isotopes (see text). 



MeV, the first derivatives, 5*271, are much closer (typi- 
cally differing by <100 keV), and, for the mixed deriva- 
tives SVpn, the differences are essentially zero. 

Figure 3 (top) shows three sets of calculations for Er 
spanning a spherical-deformed transition region. Gener- 
ally, the SkP and SLy4 functionals (both with the same 
mixed pairing) reproduce the trends in the data where 
empirical SVpn values are known (A''=84-102), including 
larger values for the spherical nuclei, a sharp drop in the 
transition region (TV ~ 90-94), and an increase beyond 
that. The lower part shows the calculated expectation 
values of /3. These are compared to empirical values 
of l/E{2+) which is a useful measure of shape corre- 
lated with the moment of inertia. The calculations repro- 
duce the transition region quite well, including the onset 
of deformation near A^=90,92 (perhaps slightly shifted 
to lower neutron numbers), and its saturation above 
iV=100. Beyond iV=102, where no data exist, the SkP 
and SLy4 results differ markedly, and in a way rather typ- 
ical for many cases. The SLy4 results often show much 
larger fluctuations than the more robust SkP results, 
which also agree better with the data than other interac- 
tions tested. Figure 3 (top) also shows SkP calculations 
where the shape was constrained to a spherical shape. 
These calculations, of course, are not suitable for transi- 
tional or deformed nuclei such as Er and their neighbors, 
but serve here as a pedagogical benchmark. Excursions 
above the spherical reference (around the middle of the 
shell) and below this line {e.g., near N=90) exhibit the 
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FIG. 3: Top: SVp„ values for Er isotopes. Calculations are 
performed with SkP and SLy4 functionals with mixed pairing. 
The SkP results constrained to a spherical shape (SkPsph) 
are marked by a dotted line. Bottom: Calculated absolute 
values of the quadrupole deformation /3 (red/green squares 
for prolate/oblate shapes) together with experimental values 
of l/£(2+). 



effects of shell structure on quadrupole correlations. In 
this context, it is interesting to recall that the microscopic 
origin of quadrupole deformations in nuclei can, within 
the nuclear DFT, be attributed to p-n interactions |23|. 

The global and local comparisons in Figs. 1-3 are com- 
plemented by a regional perspective in Fig. 4. Arguments 
[5| from generic shell structure can account for the over- 
all systematics. Normal parity orbits fill high-j, low-n 
orbits early in a major shell, but low j-high n orbits at 
the end. Thus one expects the largest SVpn values near 
the diagonal, where there is similar fractional filling of 
proton and neutron shells and therefore maximum spa- 
tial overlap. At the time of Ref. [^, there were no real- 
istic calculations that could be brought to bear on this 
empirical phenomenon. However, the power of the nu- 
clear DFT approach is shown by the present calculations 
which reproduce very well both the general magnitude of 
SVpn values and their variations across a shell, including 
larger values near the diagonal. Figure 4 illustrates an- 
other point. Overlap arguments again suggest Q larger 
empirical SVpn values where protons and neutrons are 
filling similar quadrants (lower left and upper right) and 
smaller values in dissimilar regions (upper left). It was 
speculated Q that SVpn values would also be small in 
the currently data-free lower right quadrant. The calcula- 
tions strikingly confirm this qualitative idea and enhance 
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FIG. 4: Empirical and calculated 5Vpn values for the major 
shells Z=50-82, Ai'=82-126, color coded according by magni- 
tude [red for the largest values, blue for the smallest] . (Upper 
panel taken from Ref. [^). The zig-zag lines in the theoretical 
panel enclose nuclei with known empirical SVpn values. Four 
boxes within these lines are empty since some of the nuclei 
involved are calculated to lie beyond the proton drip line and 
cannot be calculated reliably by the current DFT model. 



the need for new mass measurements in this region. 

To conclude, we have presented the first large-scale mi- 
croscopic DFT calculations of SVpn spanning the nuclear 
chart, comparing the results for different functionals, in- 
cluding different treatments of pairing, and confronting 
these calculations with empirical results. Overall, the 
agreement is impressive. This is especially significant 
since the average error of the calculated binding energies 
is several MeV. Yet, proton-neutron interaction energies 
obtained from binding energy differences generally match 
the data to well under 100 keV and, often, significantly 
better than that. The agreement is best in the deformed 
regions where the mean-field theory is expected to cap- 
ture essential physics. This level of agreement can be 
exploited to extrapolate to unknown nuclei when three 
of the four masses needed for a SVpn value are known. 
For example, using the predicted SVpn value for ^^^U, 
and the known masses for 238, 236^ g^j-^j 234rp-|^^ gives a 
predicted mass excess of -46325 keV for ^"^^Th. It would 
be interesting to test this compared to the extrapolation 



from systematics {3| of -46454 keV. This approach gives 
useful predictions for {Z,N — 2) and {Z — 2,N) nuclei 
in proton and neutron-rich regions, respectively. Finally, 
deviations of SVpn values do occur, and point to needed 
improvements in the density functionals, where specific 
effects or degrees of freedom enter. Example were noted 
at magic numbers, N=Z nuclei, and in the octupole- 
correlated Ra region. 
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